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Abstract 

The question of CP violating phases in supersymmetry and electric dipole 
moments (EDMs) is considered within the framework of supergravity grand 
unification (GUT) models with a light (~1 TeV) mass spectrum. In the mini- 



\Q , mal model, the nearness of the t-quark Landau pole automatically suppresses 

the t-quark cubic soft breaking phase at the electroweak scale. However, cur- 
f-^. ■ rent EDM data require the quadratic soft breaking phase to be small at the 

electroweak scale unless tan/3 is small (tan/? ~3), and the EDM data com- 



bined with the requirement of electroweak symmetry breaking require this 
phase to be both large and highly fine tuned at the GUT scale unless tan/3 is 



Q_i' small. Non minimal models are also examined, and generally show the same 

O ■ behavior. 

I. INTRODUCTION 

• i— I 

x 

The Standard Model (SM) of strong and electroweak interactions possesses only one 
source of CP violating phases: the phase in the CKM matrix. While the physical origin 
of this phase remains unknown, it appears to satisfactorily account for the observed CP 
violation in the K meson system, and future data from B factories and high energy acceler- 
ators will be able to test its validity for B meson phenomena. It was early on realized that 
supersymmetric (SUSY) extensions of the Standard Model allowed for an array of new CP 
violating phases, and that these phases could give large contributions to the electric dipole 
moments [EDMs] of the electron and neutron |lj, thus violating the known experimental 
bounds @||. Several resolutions to this problem have been proposed, e.g. one might as- 
sume the phases are quite small i.e. ~ 0(1CT 2 ) |],[3|], or suppress the diagrams by assuming 
the relevant SUSY particles are very heavy ||. Both these suggestions a priori appear un- 
satisfactory, i.e. the first possibility would appear to require a significant amount of fine 
tuning, and the second would imply a SUSY spectrum in the TeV domain, possibly beyond 
the reach of even the LHC. More recently, it has been realized that some cancellations can 
occur in different contributions to the EDMs in certain parts of the SUSY parameter space 
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H which has led to considerable examination of this possibility |6HT5|. In the minimal su- 
persymmetric extension of the SM, the MSSM, (which we define as the low energy SUSY 
extension of the SM obtained by supersymmetrizing the particle spectrum with two Higgs 
doublets) this does appear to reduce the amount of fine tuning needed, allowing for larger 
phases. The existence of such phases could also have effects on other predictions of the 
MSSM [16|-|21|, thus allowing future experimental checks of this idea. 

In this paper we consider these questions from the viewpoint of supergravity (SUGRA) 
grand unified models (GUTs). In particular we consider gravity mediated models with Im- 
parity invariance, where supersymmetry is broken in a hidden sector at some scale above the 
GUT scale Mq (presumably at the string or Planck scale Mpj), the breaking being transmit- 
ted to the physical sector by gravity f22[. (For previous work from this viewpoint see Refs. 
PJ?|JP,p!0|, p^ , p^1 .) Such models are much more constrained than the MSSM. Thus low en- 
ergy properties are now determined by running the renormalization group equations (RGE) 
from Mq to the electroweak scale. CP violating phases that in the MSSM are arbitrary 
parameters now get correlated by the RGE, and no longer can be assigned independently 
to satisfy the EDM experimental bounds. Further, the radiative breaking of 577(2) x U(l) 
puts additional constraints on the CP violating phases. We find that as a consequence of 
these constraints, one must still fine tune some of the CP violating phases at the electroweak 
scale to be quite small, except for small regions of the parameter space with tan/? close to its 
minimum value, or when SUSY masses are very large. Further, except for the lowest tan/3 a 
serious fine tuning develops at the GUT scale. These phenomena appear to be true both for 
the minimal mSUGRA model with universal soft breaking masses, and with nonuniversal 
soft breaking extensions. Hence the required smallness of the phases also tends to reduce the 
possibility of collateral evidence in other phenomena for their existence in SUGRA models. 

In Sec. 2, we discuss the solutions of the RGE for mSUGRA models and explain which 
phases that are large at the GUT scale get naturally suppressed at the electroweak scale 
(and which do not), and what the constraints of electroweak breaking imply both at the tree 
and 1-loop level. We then examine the amount of fine tuning needed to satisfy the existing 
EDM data, and what would be needed if those data were improved by a factor of 10 without 
finding any EDM. Sec. 3 considers some effects of nonuniversal soft breaking on the EDM 
results. 

In gravity mediated SUGRA models, the structure of the soft breaking parameters is to 
be deduced from the nature of the Kahler potential and gauge kinetic function at the GUT 
scale P2]-F21[ . Thus what might be reasonable sizes for CP violating phases are presumably 
Planck scale physics questions. While at present there is no theory of such phenomena, 
one can still examine the general framework to see under what circumstances small or large 
phases might occur "naturally" at the GUT scale. This is analyzed in Sec. 4, and a model 
with "naturally" small CP violating phases is discussed. 

Concluding remarks are given in Sec. 5. 



II. EDM FOR MSUGRA MODELS 

Supergravity GUT models with universal soft breaking of supersymmetry, mSUGRA, 
depend upon five parameters at the GUT scale: mi/ 2 (the universal gaugino mass), A 
(the cubic soft breaking mass, B (the quadratic soft breaking mass), fiQ (the Higgs mixing 
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parameter) and mo (the universal squark and slepton mass). Of these, the first four may 
be complex. However, it is always possible to make a phase transformation on the gaugino 
fields to make m.1/2 real, and since the reality is preserved by the RGE at 1-loop order, we 
will assume for now on that is real. We parameterize the remaining phases at Mq by : 

A = \A \e iaoA ; B = \B \e td0B ; /i = |// |e i6 V (1) 

The RGE allow one to determine the low energy parameters in terms of those of Eq.([l|) and 
mi/2 and m . Thus at the electroweak scale, one obtains a different A parameter for each 
fermion (A u , A d , A t , A b , A e and A T ), real gaugino masses m 3 , m 2 and mi for the SU(3)c, 
SU(2)l and U(l)y sectors, and complex B and /i parameters. 

The EDM, df for fermion f, appears in the effective Lagrangian Lf as 

Lf = --djfa^fF^ (2) 

The basic diagrams giving rise to df are shown in Fig.[l], and involve neutralino (x°, i=l,4) , 
chargino (xf, i=l,2) and gluino (g) loops with squarks (q) and sleptons (/). 

For the neutron dipole moment d n one must also take into account the gluonic operators 

L G = ——d G f a bcG a ^ a Gfr u G l l !: u (3) 



and 

c -- l d c qa^ T a qG7 (4) 
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where G^ u = ^e^^Gcap, f a b c are the SU(3) structure constants and T a = |A a , where A a 
are the SU(3) Gell-Mann matrices. In addition to the diagrams of Figjl] with 7 replaced by 
g, operators of Eqs.(|^-|4]) receive contributions from the two loop Barr-Zee diagrams of Fig|| 
|25| , and the two loop Weinberg type diagram of Fig.[3] |26| . For the neutron dipole moment 



one must use the QCD factors r] ED , r] G , i] c to evolve the results at the electroweak scale 
down to 1 GeV |27|]. 



The calculation of d n suffers from QCD uncertainties from several sources. These include: 
(1) How to relate the quark moments d u , d^ to the neutron moment d n . We use here the 



non relativistic quark model relation. (For other approaches see Ref. [|14 



1 

d n = -(Ad d -d u ) (5) 

(2) How to relate color and gluonic contributions to the electric dipole moment. We use 
here the "naive dimensional analysis" of Ref. |28|] . 

(3) Uncertainty in m s , the strange quark mass, which affects the determination of m u and 
rrid- While the quark mass ratios are well known p9f . 

Tfl Tfl 

— = 0.553 ± 0.043; — = 18.9 ± 0.8 (6) 

m d m d 

one has m s = (175 ± 25)MeV from QCD sum rules and m s = (100 ± 20 ± 10) MeV (2 GeV) 
from the quenched lattice calculation |I0| (Using unquenched approximation one expects an 
even smaller value [j3T|j ). 
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Thus the values of d n calculated below have a significant uncertainty. In spite of this, 
it is useful to see the effect on the SUSY parameter space of simultaneously imposing the 
experimental constraints on both d n and d e and we will discuss below how the uncertainty 
in m s effects these constraints. 

To specify our phase conventions, we give the mass matrices for the particles entering in 
Fig.|l]. With a convenient choice of phases, the chargino and the neutralino mass matrices 
are: 

m 2 \?2M w sinf3\ . . 

V2M w cos[3 -\n\e i6 ) [ ' 




1 mi 





a 


b 







m 2 


c 


d 




a 


c 





1 1 if 

\fi\e 


> 


V b 
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|/i|e iS 





) 



(8) 



where a = —Mzsin9wcos/3, b = Mzsin9wsinj3, c = —cot9wo>, d = —cotOwb, tan/3 = v<ijv\ 

(fl,2 



|< H\ 2 >|) and 9w is the weak mixing angle. The phase 9 is given by 

e = e 1 + e 2 + 6 lt 



(9) 



where at the electroweak scale, < >- 
matrices may be written as 



•ui 2 e iei ' 2 , and \i 



Ml 



m 



e ia "m q (\A q \ + \/2\R q e' 



-i{G+a q ) 



im q {\A q 



= |/i|e jeM . The squark mass 
\fx\R q e^ e+a ^)\ 



m 



'IB 



where the quark mass and charge, 



m 
rn 



m 2 Q + m 2 q + (I/ 2 - e q sin 2 9 w )M 2 z cos2(3 



2 

1L 

2 qR = mfj + m 2 + e q sin 2 9 w Mlcos2(3 



(10) 



(11) 
(12) 



A q = \A q \e lolq , m 2 Q and mfj are given in Ibanez et al. f32| and R q = cot(3(tan/3) for u(d) 
quarks. Similar expressions hold for slepton masses, with phases a\. Our sign conventions 



for A q and \x are those of Ref. |35 



The condition for electroweak symmetry breaking is obtained by minimizing the effective 
potential V e ff with respect to t>i, ei, v 2 and e 2 . The Higgs sector of V e ff is 



V eff = m\v\ + m\v\ + 2\Bfi\cos(9 + 9 B )v x v 2 + ^{v\ 



l2 



,2\2 



(13) 



where V\ is the one loop contribution, m 2 = fi 2 + m 2 1 . and m 2 H are the Hi 2 running masses. 
For particles of spin j one has 



1 m 2 3 

Vi = C a (-l) 2j "(2j a + l)m*(in5 - - 

64vr 2 Y J aK Q 2 2' 



(14) 



where C a is the color degree of freedom of the a th particle. In the following we choose the 
low energy scale Q to be m t (175 GeV), and include the full third generation states, t, b and 
r in V\. This allows an examination of the large tan/3 domain. 
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We view the minimization equations of V e ff as equations to determine the Higgs VEVs 
i.e. Vi, t> 2 , ei, e 2 . Thus in the tree approximation, the extrema equations dV e ff/dei = 
yield 2\B [i\sin(9 + 6b) = 0, and hence the minimum of V e ff requires 

6 = ti-6 b (15) 

(the choice of 6 = —6b leads to a maximum). The one loop corrections depend only on the 
mass eigenvalues, and hence from Eq.([L(]) only on the phase 6 + a q and 6 + ol\. Thus, at the 
one loop level, one gets a correction to Eq.flTo]) of the form [|TB].[T7 



6 = ti -6 B + fi(n -6 B + a q ,it -6 B + cn) (16) 

where f\ is the one loop correction with 6 approximated by its tree value, Eq.fllSf). As we 
will see, this correction can become significant for large tan/3. The EDMs however depend 
on both the mass eigenvalues and the rotation matrices. Hence they will depend on 6 + a q 
and a q separately, with 6 determined by the electroweak symmetry breaking condition in 
Eq.©. 

The current experimental 90% CL. upper bounds on d n |jj and d e are 

(d n )ex P < 6.3 x l(T 26 ecm; (d e ) exp < 4.3 x 10~ 27 ecm (17) 
In discussing the data, it is convenient to define the quantity 

K = log 10 \-^—\ (18) 

\flf jexp 

Thus K = corresponds to a theoretical value which saturates the current experimental 
bound, while K = —1 would represent the case, should the experimental bounds be reduced 
by a factor of 10. As pointed out in Refs. fl^Iol, various cancellations can occur among 



the different contributions to the EDM. This is illustrated in Fig.|| where K is plotted vs. 
mo for the electron dipole moment (eEDM). We note that eventually for very large mo, the 
curves fall below the K = bound (as expected). Further the allowed range of m (so that 
K < 0) decreases with increasing tan/3 (and the allowed range would become very small 
should K = —1, i.e. the experimental bounds be reduced by a factor of 10). In addition, 
the position of the dips moves to lower itlq with increasing tan/3. This happens due to the 
fact that the regions on the right of the dips are dominated by the contribution from the 
chargino diagram. As tan/3 increases, the contribution from the chargino diagram increases 
much more than the contribution from the neutralino diagram. However, a decrease in m 
increases the neutralino diagram much more than the chargino diagram. As a result the 
dips shift towards the smaller mo value for larger tan/?. We use negative values (arising from 
the factor it in the Eq. [15D of \x to satisfy the experimental constraints on the BR of b — > 57 

is. 



We begin our analysis by examining the RGE to determine what GUT scale parameters 



lead to acceptable CP violating phases at the electroweak scale. We have used Ref. |35|] for 
the RGE relating the GUT scale parameters to the parameters at the electroweak scale. In 
general, the RGE must be solved numerically and all results given below are obtained by 
accurate numerical integration. However, approximate analytic solutions can be found for 
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low and intermediate tan/3 (neglecting b and r Yukawa couplings) or for the SO(10) limit 
of very large tan/3 (neglecting the r Yukawa coupling). These analytic solutions give some 
insight into the nature of the more general numerical results. 

The solution for the A t (t) parameter in the low and intermediate tan/3 case can be cast 



in the form 5c 



A R (t) H 3 (t) 

A ° = ^ + T^) m ^ (19) 



where Ar is the residue at the t-quark Landau pole, 



and 



A R = A + m 1/2 (H 2 (t) - (20) 



D = l-Q^Y(t). (21) 



Here t = 2ln(MG/Q), the form factors E, F, if 2 arid H3 are real and are defined in [32 
and Y(t) = h t /16iT 2 , where ht is the t-quark Yukawa coupling constant. -Do vanishes at 
the t-quark Landau pole (for Q=m t , D = 1 — (m t / '200 sin j3) 2 ) and is generally small (i.e. 
Dq ~ 0.2 for m t = 175 GeV). The imaginary part of Eq.(|19|) gives: 

D \A \sinaoA = \A t \sina t . (22) 

Thus even if a$A = 7r/2, a t at the electroweak scale will be generally suppressed due to 
the smallness of Do, i.e. the RGE naturally make the phase at small due to the nearness 
of the Landau pole. (This result has been previously observed, for low tan/3 in ||.) The 
approximate SO(10) solution with large tan/3 (where Y t = Yj,) gives a similar analytic form 
for A t (t) = A b (t) with the factor 6 replaced by 7 in Eq.(pT|). Thus the suppression effect on 
a t occurs over the entire tan/3 domain. These effects thus allow «oa to be large without a 
priori violating the experimental EDM bounds. 

In contrast there is no analogous RGE suppression effect in the first generation A pa- 
rameters A u , Ad and A e . From the RGE, one finds that A u has the form of Eq.fliDD for low 
and intermediate tan/3 with the factor 6 replaced by 1 in Eq. (pip . Now D Q ~ 1 and no sup- 
pression effect occurs for a u at the electroweak scale in the analogue of Eq. . However, 
the EDMs are not very sensitive to the first generation a u , a^, a e and one can have large 
values of these parameters without violating the experimental bounds in a reasonable region 
of the SUSY parameter space. 

From the RGE, one sees that the phase of the fi parameter does not run at the 1-loop 
level, i.e. = 9^, where 9^ is the phase at the electroweak scale. Thus a large GUT scale 
phase will lead to a large electroweak scale phase. However, as seen in Eqs.(|7|-|roD, 9^ enters 
only in the combination 9, and 9 is determined in terms of 9b, ai and a q by the minimization 



of the effective potential, as given in Eq.(16). Thus (with our choice of the phases in the 
mass matrices) 9^ does not enter in the EDMs. 

The EDMs, however are highly sensitive to 9b, which by Eq. (|16D , enters into all the 
diagrams, and over most of the parameter space, 9b must be small to satisfy the EDM 
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experimental bounds. This is illustrated in Fig.|5] where we have plotted the eEDM K vs. 
9 B . We see that 9 B can be large i.e. 9 B ~ 0.08, for low tan/? = 3, though it becomes smaller 
for higher tan/?. However, even for tan/3 = 3, the range A9 B of 9 B so that K < is very 
small i.e. A9 B ~ 0.02 (and 9 B would become quite fine tuned if the experimental bound 
were reduced by a factor of 3 i.e. K = —0.5). Figs. 6 show that A9 B remains small for 
mi/2 ~ 350 GeV (rrig ~ 1 TeV) and only relaxes somewhat when mi/ 2 — 700 GeV [rrig ~ 2 
TeV). In this figure we have plotted the K < regions as a function of 9 B and m\/2 for 
three different values of tan/3. Since the contribution from the chargino diagram increases 
with tan/3, the larger tan/3 requires a smaller phase to produce the necessary cancellation 
between the chargino and the neutralino diagrams. 

The RGE allow us to examine the significance of the above results at the GUT scale. 
The RGE for B and A t can be solved in the low and intermediate tan/3 regime to obtain 

B = B - i(l - D )A - $m 1/2 (23) 

where 

* = ~(1-A>)fr + [3^ + ^i]^f (24) 

hi = t/(l + Pit), «g = 1/24 is the GUT gauge coupling constant and $ are the MSSM beta 
functions. The real and the imaginary parts of Eq . (^3|) give: 

\B\sin9 B = \Bo\sin9oB — tt(1 — D )\A Q \sinaoA (25) 



\B\cos9 B = \B \cos9 0B - -(1 - D )\A \cosa 0A - $mi/ 2 (26) 
Eqs.(||j2|) can be viewed as determining the electroweak scale values of \B\ and 9 B in 



terms of the GUT input parameters. Thus 9 B depends upon both the initial phase 9q B 
and the A Q phase a A- Alternatively, one may use Eqs. (|25| , |26|) to impose electroweak scale 



phenomenological constraints on the allowed GUT parameters. The requirement that the 
GUT theory gives rise to electroweak symmetry breaking gives 

\B\ = ^in2(3-^ (27) 

where ml = a\ + /i^ and = |/i| 2 + m 2 H . + Ej. (Here m 2 H . are the running Higgs masses 
at the electroweak scale, \fi\ is determined by electroweak breaking and Ej are the one loop 
corrections [37].) In addition, we have seen that the experimental bounds on the EDMs 
restrict 9 B to be small, i.e. \9 B \ ~ 0.1 (and usually much smaller) and the allowed range, 
A9 B , that satisfies the EDMs is quite small. These conditions put severe constraints on the 
GUT scale theory which we now discuss. 

Consider first the case where tan/3 is very close to its minimum value e.g. tan/3=2. 
Here l/2sm2/3 is not small and by Eq. (|27|) , \B\ is of normal size. Eq.(p5|) then implies 



that 9 B and 9q B are of roughly the same size (even for o>qa — tt/2) as previously noted 
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in Ref. ||. Thus a reasonable GUT theory requires one only to justify the size of 9qb- 
However, as tan/3 increases, l/2sin2(3 decreases, and unless /i becomes anomalously small 
or m| anomalously large, Eq.fl2T|) implies that \B\ becomes small. Thus the lefthand side of 
Eq . (|25"D becomes small (and to first approximation can be neglected), and Eq.(^) implies 
that 9 b is determined by ctoA mainly, i.e. 9ob is large if o>oa is large. This result is illustrated 
in Fig.[7| where 9 b is plotted vs. 9b- We see that since 9b is small for all tan/3, 9 b is in 
general large, e.g. even for tan/3 = 3 one has 9qb — 0.8. 

Returning to Eq. (|25|) , we see, however, that A9 0B , the range of values of 9 b that will 
satisfy the EDM constraints, is very small. Thus for fixed A , Eq.fl2"5]) implies for large tan/3 
that 

A9 0B = |tA0b < A9 B (28) 
I -^0 1 

and since A9b is small, A9qb will be very small. This is illustrated in Fig.[8] where the 
allowed values of 9 b for K < are plotted vs. for d e and d n for tan/3=3, 10, 20. 

We see that the phenomenological constraints at the electroweak scale imply that 9qb is 
both large and its value is sharply fine tuned (unless aoA is small, SUSY masses are large 
or tan/3 is small). One may alternatively view this from the "top down": if Planck physics 
determines Oqa and 9 b to be large and fixed, then if the model is to achieve electroweak 
symmetry breaking with EDMs below existing bounds, there will be a fine tuning of other 
GUT scale parameters (unless tan/3 is small or the SUSY masses are large). 

As discussed above, there are a number of uncertainties in the calculation of the neutron 
EDM. Fig.|9| illustrates the effects of varying the d-quark mass in a plot of the region allowed 
by the experimental EDM bounds (K < 0) in the mo — m.1/2 plane for 9b = 0, \Aq\ = 300 
GeV, aoA = tt/2, tan/3=3. The bound on d e already excludes all regions below the lower 
solid curve, while the d n bound excludes regions to the left of the upward running curves for 
m d = 5 MeV (dotted), m d = 8 MeV (solid), m d = 12 MeV (dashed). Thus the combined 
exclusion region increases from ~ 260 GeV to m 1( / 2 ~ 440 GeV as m d increases, a 
significant change. In the other figures in this paper we have used the middle value of 
rrid = 8 MeV corresponding to m s ~ 150 MeV. 

In Fig.p, we have chosen 9b = 0. As one increases 9b, one will eventually arrive at a 
neutralino-chargino cancellation region, and this will reduce the values of mo and m\i% which 
are excluded. Further, for low tan/3 this cancellation can occur at relatively large 9b (i-e. 
9b — 0.2) as shown in Ref. 0. Alternatively, increasing tan/3 can also cause this cancellation, 
as the loop corrections, Eq . (|i~6D , become large with large tan/3 and the contribution from the 
chargino diagram increases with tan/3. This is illustrated for tan/3=20 in Fig {Tl] where the d n 
curve for K = (solid) bends downward below m =500 GeV when 200 GeV< < 400 
GeV, showing in this exceptional type situation that one can satisfy the EDM constraints 
with a light particle spectrum and large tan/3. Note that the cancellation can occur for a 
relatively wide band of m^, and persists even for K = —0.5. 

The regions of the SUSY parameter space that get eliminated by a joint imposition of 
the experimental bounds on the d e and d n are sensitive, of course, to the choice of the 
parameters. Thus in Figs.|6| one sees that for m = 100 GeV, one requires ~ 350 

GeV (nig ~ ITeV) to jointly satisfy the d e and d n bounds, and this requirement is roughly 



independent on tan/3. Fig. 11 shows that if is increased to 250 GeV one requires now 



8 



mi/2 ~ 160 GeV (m§ ~ 450 GeV) to satisfy both bounds, i.e. if m is increased, the lower 
bound on m 1( / 2 is decreased as expected. However, if \Aq\ is increased, it is not necessary that 
the allowed domain of mi/ 2 increases. This is illustrated in Fig.|12] for the choice |y4 |=800 
GeV, m =100 GeV. The increase of \A Q \ raises the allowed values of 9 B for both d e and d n , 
but the raising of d e is larger causing the joint allowed domain to require m ly / 2 ~ 550 GeV 
(mg ~ 1.5 TeV). 

In order to exhibit the size of | fi | , we have plotted | function of mi/ 2 for m =100, 

300 and 700 GeV in Fig.|T5|. The EDM constraints have not been imposed in Fig.|T3], but 
from Fig. 6a, one sees that they are satisfied at least for m =100 GeV for d e for essentially 
the entire mi/ 2 range (which allows 9 B ~ 0.08). Note that | // | is large (i.e. | fi | 2 /M| 1) 
over the full domain so that one is generally in the gaugino scaling domain. 



III. NONUNIVERSAL SOFT BREAKING 



There are three types of nonuniversalities that might be considered in SUGRA models: 
nonuniversal gaugino masses at Mg, nonuniversal scalar masses at Mg and generation off 
diagonal cubic soft breaking matrices Ai~ and scalar mass matrices mf ■ . We will discuss here 



briefly the first two possibilities, a more detailed analysis will be given in a later paper [38 
(The third type of non universality has the possibility of generating e and e'/e CP violations 
independent of the CKM phase [p9"| .) 

In general, the U(l), SU(2) and SU(3) gaugino masses at Mg can have the form 



m 1/2 i 



|«ii/ 2 ;|e^ 1 ; i 



1,2,3 



(29) 



A convenient phase choice is to set the SU(2) phase to zero: <p 2 = 0, and hence the chargino 
mass matrix, Eq.(|7]), is unchanged. The remaining phases produce effects in the RGE for 
Ml): 



A t = DqAq + ®ti\mx/2i\e 



i<j>i 



(30) 



and hence the imaginary part gives 

\A t \sina t = \A \D sina 0A + $ 4i |m 1/2 i|sm0i 



(31) 



Thus the Landau pole factor D still suppresses any large phase aoA at the GUT scale, 
though a t may become large due to the <pi phase which does not run at one loop order and 
hence is not suppressed. This phenomena can affect d n . 

A second effect occurs in the neutralino mass matrix which now reads 



M x o 



( 


mi e* 





a 


b \ 







rh 2 


c 


d 




a 


c 





|/i|e* e 


V 


b 


d 


|/i|e ie 


/ 



(32) 



where 9 = e\ + e 2 + 0«. Thus (f>i effects any neutralino-chargino cancellations in d e and 



d n . An example is illustrated in Fig. 14. One sees that 9b can be large but is still tightly 



9 



constrained to satisfy the experimental EDM bound (i.e. A9b — 0.02). One is again led to 
a fine tuning at the GUT scale. 

As an example of nonuniversal sfermion masses, we consider an SU(5)-type model where 
the GUT group possesses an SU(5) subgroup with matter embedded in the usual way in 
10 + 5 representations. We also assume the first two generations remain universal to suppress 
flavor changing neutral currents. The Higgs masses have the form at Mq 

m ik = m l{ 1 + <*i); m 2 H2 = ml(l + S 2 ) (33) 

while the third generation 10 representation (containing ti, £r,6z,, ?r) and the 5 (containing 
bn, v T -> ?l) masses are parametrized by 

m 2 10 = m 2 (l + 6 10 ); = m 2 (l + 6 B ) (34) 

Here the 8i > — 1 represent the deviations from the universality. In addition the third 
generation cubic soft breaking parameters are A ot and A ob = A 0t . 

One may use the non universal sfermion and Higgs masses to soften the effects of the 
experimental EDMs. This is illustrated in Fig.|l5| where the K=0 constraint is imposed in 
the m — mi/2 plane for a choice of Si parameters. For Si = 0, the d n curve would continue 
to rise for small mjw as seen in Fig.^, rather than turn over as in Fig.[0]. For example, 
for m 1 / 2 =220 GeV and tan/5=5, K(d n ) < occurs for m > 580 GeV in the non universal 
case, whereas K(d n ) < occurs for m > 750 GeV in the universal case. The particular set 
of SiS in the Fig.[T5] is chosen to reduce the stop and the sbottom mass and to satisfy the 
requirement of the radiative electroweak symmetry breaking. The lower stop and sbottom 
mass increases the Weinberg type diagram which has a subtractive effect in the net nEDM 
magnitude. 



IV. SUGRA MODEL OF SMALL PHASES 

While choosing the CP violating phases to be small, i.e. (pi ~ 10 -2 , to satisfy the EDM 
constraints appears to be artificial in the low energy MSSM, it is possible that a natural 
model of this type can arise in SUGRA GUT models. 

SUGRA models depend on three functions of the chiral fields $ Q : /«/?($«) the gauge 
kinetic energy, K(§ a , ) the Kahler potential and W($ a ) the superpotential. One assumes 
a hidden sector exists where some fields, e.g. moduli or dilaton, $j grow VEVs of Planck 
size to break supersymmetry i.e. 

Xi = k < $i >= 0(1) (35) 

where k~ 1 = Mpi = 2.4 ■ 10 18 GeV. Thus one can write = $ a } where <3? a are 

the physical sector fields. One might expand the Kahler potential in a power series of the 
physical fields: 

K = K- 2 C^ + (C^^ b + ic2,*a*k*c + ...) + (eg*,,*! + ^M^c + -) (36) 

where M is a large mass. The are dimensionless and we assume them to be 0(1). The 

first parenthesis is holomorphic, and can be transferred to the superpotential by a Kahler 
transformation giving rise to the \x term: 
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W _ W + «a W ( c W$ a $ 6 + — c2i$ a $ 6 $ c + ...) (37) 

The leading additional terms on the right occur when W there is replaced by its VEV 
(k 2 < W >~ 777.3/2) and one of the fields in the cubic term, e.g. $ c , has a GUT scale VEV. 
Then W -> W + /j, ab & a &b, where 

Vab = + ^C$ c )m 3/2 (38) 

In the perturbative heterotic string we expect M ~ Mpi and hence Mq/M ~ O(10~ 2 ). If 
one assumes is real and c^ 3 - 1 is complex with arbitrary size phase, then jx naturally has 
a phase of size 9^ ~ 0(1O -2 ). A similar analysis can be done for the other parameters (A, 
B, to) yielding automatically small phases. One assumes only that the leading (renormaliz- 
able) term in the matter expansion is real, while the higher terms (presumably arising from 
integrating out the tower of massive states) can have arbitrary phases. 

Other string scenarios can have M ~ 0(Mg). In that case one expects the nonrenor- 
malizable terms to give rise to large phases, unless the low lying members of the tower of 
heavy states do not couple to the fields of the physical sector. 



V. CONCLUSIONS 

In SUGRA GUT models, the basic parameters are given at the GUT scale. The renor- 
malization group equations then carry this information to the electroweak scale, allowing 
for experimental tests of the model. Thus what represents a "natural" choice of parameters 
is presumably a GUT scale question. 

The experimental values of the electric dipole moments put strong constraints on the low 
energy SUSY parameter space. However, for mSUGRA models at least, the RGE suppress 
the value of the A t phase at the electroweak scale, at, so that any size of phase, a^A, even 
7r/2 at Mg will generally lead to acceptable EDMs. This effect is due to the nearness of the 
t-quark to its Landau pole. 

The EDMs, however, are very sensitive to the B phase at the electroweak scale, 9 B . In 
the domain of light SUSY mass spectrum (i.e. ~1 TeV) and a A — 0(1), the EDMs allow a 
large 9b (i.e. O(10 -1 )) only for low tan/5, i.e. tan/5 ~ 3. Further, for fixed SUSY parameters 
in this domain, the allowed range A9b of 9b needed to satisfy the EDM constraints is very 
small for tan/5 ~ 3. 

The above results strongly affect the GUT theory. Thus for a A large and a light mass 
spectrum one finds at the GUT scale that 9 0B is generally large, i.e. 0(1) even for low tan/5, 
and the allowed range A9 0B for satisfying the EDMs is exceedingly small for tan/5 ~ 3. Thus 
for this situation, 9qb must be chosen large and be very precisely determined as tan/5 gets 
large, leading to a new fine tuning problem at the GUT scale. The origin of this difficulty 
resides in the requirement that the GUT models give rise to parameters at the electroweak 
scale that simultaneously satisfy the experimental EDM constraints and give rise to radiative 
electroweak breaking. 

The above discussion suggests that reasonable mSUGRA models with large phases and 
mass spectrum below 1 TeV can only be constructed for small tan/5 i.e. tan/5 ~ 3. The 
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current LEP data combined with cosmological constraints already require tan/3 rC 2 pi, and 
LEP, the Tevatron and the LHC will be able to explore the domain tan/5 ~ 50 ^T|. Thus 
the value of tan/3 (or a lower bound) is something that will be determined experimentally in 
the not too distant future. Should tan/5 turn out to be large, it is possible that only models 
with small GUT phases 4> 0i ~ O(10~ 2 ) are reasonable, and it was shown above that such 
models can arise naturally with O i — 0{Mq/Mp{). 

Our analysis has included the loop corrections to the radiative breaking condition. While 
these corrections are small, they grow with tan/? and are competitive to the 9b contribution 
to the EDMs. Thus for large tan/3 (e.g. tan/5 ~ 20) they can produce a new cancellation 
phenomena (with #b=0), allowing the EDM constraints to be satisfied for a light mass 
spectra (e.g. m ~300 GeV, m 1 / 2 — 200 GeV) which for lower tan/5 would be excluded by 
the EDM data. 

The combined restrictions of the experimental EDM bounds for the neutron and elec- 
tron put strong constraints on the SUSY parameter space. However, since the theoretical 
cancellations needed to satisfy experiment are delicate, as one varies the parameters one 
can get rather different excluded regions. Thus, while we have found that increasing m 
gives an allowed region of smaller m\/2 (as one might have expected), increasing | Aq | can 
require m.1/2 to be increased to satisfy the EDM bounds. However, it should be stressed 
that there are significant theoretical uncertainties in the calculation of d n , and we have seen, 
for example, that the current uncertainty in the quark masses m u and (arising from the 
uncertainty in m s ) can lead to a significant variation in the domain the neutron EDM data 
can exclude. 

Much of the above discussion holds also for non minimal models with non universal 
gaugino and sfermion soft breaking masses. A more detailed discussion of the effects of non 
universalities will be given elsewhere 
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FIGURES 




f gaugino f 

FIG. 1. One loop diagram. The photon line can be attached to any charged particle. 
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FIG. 2. Two loop Barr-Zee diagrams. A is the CP odd Higgs and qi are the mass diagonal 
squark states. 
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FIG. 4. K (defined in the text) is plotted as a function of mo for d e . The solid, dashed and 
dotted lines are for tan/3 = 3, 10 and 20 respectively. The other input parameters are aoA = § , 
\Aq\ = 300 GeV, 6 B =0.02 and m 1/2 = 300 GeV. 
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FIG. 5. K (for eEDM) is plotted as a function of 9b- The solid, dashed and dotted lines are 
for tan/3 = 3, 10 and 20 respectively. The other input parameters are uoa = f , |A)| = 300 GeV, 
m 1/2 = 300 GeV and m = 100 GeV. 
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FIG. 6. The K=0 contours plotted as a function of 0b and m\/2- The solid lines are for d n 
and the dotted lines are for d e . 
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a) m =100 GeV, |A |=300 GeV, tanp=3, a 0A =7i/2 
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b) m =100 GeV, |A |=300 GeV, tanp=10, a 0A =7t/2 
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c) m =100 GeV, |A |=300 GeV, tanp=20, a 0A =7t/2 
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FIG. 8. Oqb vs mi/2- Upper and lower lines are the allowed range so that K < 0. The solid 



lines are for d n and the dotted lines are for d e . 
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FIG. 9. The K=0 contours are plotted as a function of mo and m^j. The dotted, solid and 
dashed lines are for (1 GeV)=5, 8 and 12 MeV respectively. The other input parameters are 
«0A = f , | A) I = 300 GeV 6b=0 and tan/3=3. Excluded regions are below d e and to the left of d n 
curves. 
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FIG. 10. Allowed region in tuq — m 1 / 2 plane for d e and d n . The other input parameters are 
a 0A = §> | A) | = 300 GeV, 9b=0 and tan/3=20. The solid lines are for K=0 and the dotted lines 
are for K=-0.5. 



25 




26 



.25 



0.2 



. 15 



. 1 



. 05 




400 



500 



600 



700 



800 



m 1/2 (GeV) 

FIG. 12. The K=0 contours are plotted as a function of 0b and mi/ 2 , where «oa = f 5 \Aq\ 
GeV, mo=100 GeV and tan/J=3. The solid lines are for d n and the dotted lines are for d e . 



800 



1200 




200 



300 



400 



500 



600 



700 



800 



m 1/2 (GeV) 

FIG. 13. The parameter \fi\ is plotted as a function of m\j2- The three solid lines are for 
mo = 700 (top), 300 (middle) and 100 (bottom) GeV. The other input parameters are oqa = f , 
\A \ = 300 GeV, tan/?=3 and 6 B = 0.08. 
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FIG. 14. K (for d e ) as a function of 6b- The other input parameters are mo=200 GeV; 
m 1/2 j=150 GeV, tan/?=3, a A = f , |A)| = 100 GeV and a = 4vr/5 = 3 . 
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FIG. 15. Allowed region in mo — m^ 2 plane (K < 0) in a nonuniversal scenario. The other 
input parameters are oqa = ^, 6b = 0, \Aq\ = 300 GeV, tan/3=5, 5\=1, <52=-l, 5\q=-1 and 55=-0.3. 
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